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Abstract. 

The last decade has seen increasing efforts to circumscribe and bound the 
cosmological Hubble flow in terms of model-independent constraints on the 
cosmological fluid — such as, for instance, the classical energy conditions of general 
relativity. Quite a bit can certainly be said in this regard, but much more refined 
bounds can be obtained by placing more precise constraints (either theoretical or 
observational) on the cosmological fluid. In particular, the use of the iu-paramctcr 
(w = p/p) has become increasingly common as a surrogate for trying to say something 
about the cosmological equation of state. Herein we explore the extent to which 
a constraint on the w-parameter leads to useful and nontrivial constraints on the 
Hubble flow, in terms of constraints on density p(z), Hubble parameter H(z), density 
parameter £l(z), cosmological distances d{z), and lookback time T(z). In contrast 
to other partial results in the literature, we carry out the computations for arbitrary 
values of the space curvature k G [—1, 0, +1], cquivalcntly for arbitrary Slo $ 1- 
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1. Introduction 

The classical energy conditions of general relativity [U [2], despite their well-known 
limitations [31 El [5], are nevertheless very useful surrogates for controlling the extent 
to which one is willing to countenance extreme and unusual physics. As applied to 
cosmology, in addition to the very general cosmological singularity theorem presented 
in [1], the classical energy conditions have (at the cost of additional hypotheses) 
been used to place more precise limits on the expansion of an idealized FLRW 
universe [61 [7J [H El [10]. More recently, these ideas have been extended in various 
ways in [HI [121 EH EH]- In all of these analyses there is a trade-off between the precision 
and generality of the constraints one obtains — the art lies in choosing a form of the 
input assumptions that is as general as possible, but not too general, for the precision 
of the output constraints one wishes to derive. 

In the current article we shall derive some very general bounds in terms of 
assumptions about the w-parameter, where as usual w = p/p. Specifically, we shall 
ask the question: If we know for theoretical reasons, or can observationally determine, 
that w lies in some restricted range 

w(z) G [w-,w+], (1) 

between redshift zero and redshift z, what constraint does that place on the cosmological 
expansion? We shall see that considerable useful information can be extracted regarding 
the density p(z), Hubble parameter H(z), density parameter fl(z), various cosmological 
distances dx{z), and lookback time T(z). Specifically, for some generic cosmological 
parameter X(z), we shall be looking for bounds of the form 

X(z)^X f(n ,z). (2) 

Conversely, observational constraints on these cosmological parameters can be used 
to infer features of the cosmological fluid in a largely model-independent manner. In 
contrast to other partial results scattered throughout the literature, we carry out the 
computations for arbitrary values of the space curvature k G [— 1,0, +1], equivalently 
for arbitrary Qq ^ 1. 

2. Strategy 



Our strategy will be to adopt a standard FLRW cosmology 

dr 2 
1 — kr 2 

then, (setting 8tiGn — > 1, but explicitly retaining the speed of light c), we have the two 
Friedmann equations: 



ds 2 = -c 2 dt 2 + a{t) 2 { - + r 2 (d6 2 + sin 2 9 d0 2 ) \ , (3) 



a 2 kc 2 



a 2 kc 2 
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Together, these two equations imply the standard conservation law: 

P = -3(p + p)^. (5) 
We also have the fundamental definitions 



PHubble 

and 



a 2 



a 2 



3H 2 , (6) 



= _±_ _p_ R2 + k ° 2 / a2 = i , J^L (7] 

PHubbic 3H 2 H 2 + a 2 H 2 ■ { ] 

For intermediate steps of the calculation we shall work with the very simple linear 
equation of state 

p = w*p, (8) 

where u>* is taken to be a constant. Picking some generic cosmological parameter X(z), 
we shall first calculate X Wt (z), and then (by assuming that w(z) G [w-, w + ] from redshift 
zero out to redshift z, and depending on the direction of the relevant inequality) use 
this to derive bounds of the form 

X w _(z) <X(z) <X w+ (z), (9) 

or 

X w+ {z)<X{z)<X w _(z). (10) 

We shall also make the extremely mild assumption that the density is positive 

p>0. (11) 

This is certainly a completely redundant assumption for k = and k — +1 FLRW 
universes. Only for k — — 1 universes does this provide the extremely mild additional 
constraint H > c/a, that is, H(z) > (c/a ) (1 + z). [§] 

3. Density 

We now apply this strategy to the density. From 



we have 



p = -3(p + p)~ = -3p(l + w*)~, (12) 
a (X 



£ = -3(1 +«;,)-. (13) 
p a 



| Historically it was common to refer to this quantity as the "critical density", ^critical, but with the 
advent of widespread acceptance of a nonzero cosmological constant, or more generally "dark energy" , 
the logical connection between this "critical" density and possible re-collapse of the universe has been 
severed. In a modern context then, it is inappropriate to refer to this as a "critical" density, and the 
considerably more neutral phrase "Hubble density" is preferable. 

§ This is equivalent to enforcing d > c, for a k = — 1 FLRW universe, noting that d is not a physical 
velocity, so that d > c is a perfectly acceptable physical statement. 
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So integrating, for constant we obtain the well-known result 

p w , = P o(a/a o r 3 ^ = po (1 + zf( 1+w *\ (14) 

But now ask what happens if we only know that u>_ < w(z) < w + l (Where in the 
real observable universe w(z) certainly need not be a constant.) Following the above 
analysis, we find that we must replace equalities by inequalities and so deduce 

po (1 + z)^ 1+w -^ < p(z) < po (1 + z) 3(1+lu+) . (z > 0). (15) 

Note that for z > we are looking into the past; in contrast for —1 < z < 0, we are 
looking into the future [15] . and the inequality reverses to[|J 

po (1 + z) 3(1+w+) < P(z) < po (1 + zf 1+w ^; (-K z < 0). (16) 

Of course, these simple constraints on the density are by far the most elementary of the 
inequalities we shall deduce — some of the other inequalities derived below will prove 
to be much more subtle. 

If we now in addition relax our initial constraint on po, by assuming we only know 
that the present epoch density lies in some bounded interval 

Po e [p _,Po+], that is, p _ < Po < Po+, (17) 
then these two bounds generalize to 

p _(l + ^) 3(1+u '- ) < p(z) < p 0+ (l + z)^ 1+w ^; (z>0); (18) 

p _(l + < p{z) < p 0+ (l + ^) 3(1+w - ) ; (-Kz< 0). (19) 

4. Density parameter 

We have the following identity 

a 2 H 2 clqHq a 2 H 2 ll p 

This leads to the useful result 

n(z) \ n J p( z y 1 ) 

Therefore, a bound on p(z) automatically implies a bound on Q(z). From the result for 
Pw*(z) presented above, we deduce that bounds on Q(z) can be given in terms of 

- ( q + z )-V*+i) t (22) 



which we can equivalently recast as 



n (i + z) 



Q w (z) = "" uv ^ ' ~' o — tt. (23) 

We can now use this quantity, which was derived for strictly constant w*, to bound the 
density parameter Q(z) for a more realistic matter model satisfying the milder condition 
W- < w(z) < w + . We obtain: 

|| Furthermore, note that there is no reason to ever go below z = — 1, as z = — 1 corresponds to infinite 
expansion. Also, note that the sign of 1 + iu_ and 1 + w + does not affect these inequalities. 
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• If f2o < 1 (but remember that by assumption Q > 0) then 

n w _(z) <n(z) <n w+ (z); (z>Q), (24) 

V w+ {z) < tt{z) < n w _(z); (-l<z<0). (25) 

• If Q = 1 then \/z : fi(z) = 1. 

• if n > 1, 

ft w+ {z) <Q{z) <Q w _{z); (z>0), (26) 

Q w _{z) < Q{z) < Q w+ (z); (-1<2<0), (27) 

but note that the bound can break down when the denominator of Q Wr (z) equals 
zero — this occurs at 

Qo _ lX 1/(3^+1) 



znK,O )=(^— 1 -1. (28) 

The failure of the bound might occur either in the past or the future depending on 
the value of w;*. 

— If 3w* + 1 > then the bound is useful only for z > zq(w*, Qq) < 0, implying 
a limitation in the past. 

— If 3m;* + 1 = then the bound is valid for all z. 

— If 3m;* + 1 < then the bound is useful only for z < zq(w*, Qo) > 0, implying 
a limitation in the future. 

Note that nothing unusual need happen to the universe itself at zqIw*, Qq), it is 
only the bound that loses its predictive usefulness. Combining these observations 
we see that for tt > 1 it is better (in the sense of reducing the amount of special 
case exceptions to the general rule) to recast the bounds in the form: 

$r) f 1 + *>- (3 " ++1) * ^ * ( V) (1 + z) " ,3 "" + " (29) 

for z > 0, and 

(1 + ^- +,, * %r £ ( V) (1 + zr<3 " ++ " (30) 

for -1 < z < 0. 
5. Hubble parameter 

Let us now use the density equation (the first Friedmann equation) and the definition 
of the density parameter Q to write 



h> = £ - k 4 = > * - 4^ = > tufl - 4(n - ijflj. pi) 

3 er p 3 or Oq po a 
That is, as an identity: 



H 2 = Hi \ Q -P- - (O - 1)4) = #0 ffio- " ("o - 1)(1 + zf \ . (32) 
Po a 2 ) { po 
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But we have already derived a formula for p Wt (z) , whence 

HlSz) = Hi {no(l + zf^ - (fi - 1)(1 + zf) , (33) 
which we can recast as 

H Wm (z) = H (l + z )^/l-n + n (l + zf w *+\ (34) 
For realistic matter, satisfying some constraint u>„ < w(z) < w + , we then deduce 

H w _(z)<H(z)<H w+ (z); {z > 0); (35) 

H w+ (z)<H(z)<H w _(z); (-K z < 0). (36) 

Note that the Hubble bound ceases to provide useful information once the argument of 
the square root occurring in H Wt (z) becomes negative. 

• For Q < 1 there is no limitation in the physical region z G (—1, oo). 

• For Qo > 1 this limitation manifests itself at Zh(w*,Qo) = zq(w*, f2 ), the same 
place that the bound on Q(z) ran into difficulties. (Some numerical estimates of 
where the bounds fail, based on current consensus observational data, are discussed 
in 



Finally, suppose that we do not have precise information regarding H and Qo, and only 
have the more limited information 

H e[H _,H 0+ ], n e[n _,n 0+ ], (37) 

then these two Hubble bounds further generalize to 



flo_(l + z)Jl-fi _ +a>_(l + ^) 3 ™- +1 < H{z) 



< H 0+ (l + z)Jl - Q 0+ + fio + (l + ^) 3 ™+ +1 ; (z > 0); (38) 



H _(l +z)Jl- n 0+ + Q J1 + z)^ 1 < H(z) 



< H 0+ (l + z)yjl -Vl _ +n _(l + zf w -+ 1 ; (-Kz<0); (39) 
subject to the caveat that for Q > 1 we should not push the bound past zh(w*, fl ). 

6. Cosmological distances 

Let us for the time being focus on Peebles' definition of "angular diameter distance". 
This is what Weinberg calls the "proper motion distance" [HI HE], for more definitions 
and a discussion regarding the physical interpretation of the cosmological distance scales 
see [H] , see also [H3, 02] . We make this choice to minimize the number of factors of 
1 + z in subsequent formulae. Then the standard definition is 

dp = a sin fc ( — °— [ -^rdz) . (40) 



aoH J H{z 
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a H 

this can be rewritten more suggestively as 
c 1 



sinh 



o 



1 — On 



Hp 
H(z) 



dz 



(41) 



(42) 



When interpreting this last formula for Q > 1 we make use of the fact that sinh(iO) = 
zsin(O). Substituting H(z) — > H Wt (z) and performing the integral, after considerable 
effort both Mathematica and Maple yield 



Hp 

H Wt {z 

whence 

dp,., (z) = 



-dz = In . 



, (43) 



-2/(3w,+l)' 



(yr^r^ + i) (1 + 2 



,(3™*+l)/2 



v^T^+ v/l-fio + fio(i + 2)( 3 ^ +1 ) 

This simplifies slightly 



(44) 



dp... (z) 



2Hq\/1—TIq 
- (1 + z)- 1 



[1 + z) 



(vT^+i; 



2/(3iu,+l) 



v / T^T^+ v / l-^o + fio(l + 2)( 3 ^+ 1 ) 

, \ -2/(3m,+l)' 

(VT^+1) 



We now note 



(45) 



,(46) 



VT^TTo + v / l-^o + fio(l + ^) (3u '* +1) (1 - VT^) (1 + 

(cross multiply top and bottom), which finally permits us to write the most tractable 
form of our exact result for Peebles' angular diameter distance (in a constant w(z) = 
FLRW universe): 

2/(3iu*+l) 



dp., (z) 



2H ^T=W Q (1 + z) 



1 - fto + fto(l + z)^ w *+ l ) - VT^TTo 



2/(3iu*+l)' 



(47) 



In this final expression we are always raising quantities to the same power, and the 
difference between the two terms is just in the placement of + and — signs. (Note that 
this expression is guaranteed to be real whatever the value of Qq; for Q > 1 the two 
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terms are complex conjugates of each other and after taking the pre-factor a/1 — f2o hito 
account, the overall combination is guaranteed to be real.) 

Note that once we have an explicit formula for the (Peebles) angular diameter 
distance dp, any of the other standard cosmological distances can easily be obtained by 
multiplying by suitable powers of (1 + z) [UJ [TBI HE], see also [13 EE]. In particular the 
luminosity distance is 

2/(3tu*+l) 



dL Wt (z) 




1 - fi + fto(l + z) (3w * +1) + 



2/(3w,+l)' 



l + v^T^) 



(48) 



Returning to Peebles' angular diameter distance, the Taylor series expansion in z can 
be computed as 



d P w A z ) = Jf Q \ z 



Q + 3w*Q 2 
- z 



+ 



+ w*(2fi 



en§) 



-6tt + 9Q 2 ) 



z 3 + 0{z 4 



(49) 



Perhaps of more interest is the Taylor series expansion in Q (since observationally we 
have good reasons for expecting Q « 1). The leading term is easy to calculate 

2c T {l - (i + zy^+v} + o[n - 1]. 



dp™, 0) 



i?o(3tf* 

Extracting the next 0[Qq — 1] term is not too difficult, but is somewhat tedious 

d P M= rrj\ U {!-(! + *)- 

H (3w* + 1) 



(50) 



-(3uj.+1)/2 



} 



[go ~ 1] C 

Ho 
1 

" 6 



- (1 + ^)-(3^+l)/2 

(3w* + 1) 

(1 + z )"(3«U+l)/2 



,3(3to»+l)/2 



3(3w* + 1) 



(3w, + l)/2 

+ O([fi -1] 2 ). (51) 

In any realistic situation (provided you accept the standard consensus cosmology) the 
uncertainties in w will completely dwarf any possible effect due to uncertainties in fi , 
so carrying the expansion to higher order is not warranted. 

As usual, d Pwt (z) [or dL Wt (z)] can be used to bound dp(z) [or d^z)]. Specifically, 
let w lie in the range [«?_, w + ] then independent of f2 ^ 1: 

• <Zp„, + (*) <dp(*) <d Pw _(z); (z>0), 

• dp„_(z) < dp(z) < d Pw+ (z); {-Kz< 0). 
Here dp w (z) is given by the rather formidable equation (T47J). 
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Finally, consider the "lookback time" defined by: 



T(z) = 
That is: 



dt 



dt = / —da = / - — 



da 



a da 



a a 



1 d(ao/(l + z)) 
H a /(l + z) 



1 dz/(l + z) 
H 1/(1 + z) 



T(z) 



dz. 



(1 + z) H(z) 
Using the known form of H Wt (z) we define 

T„. (z) = — 



1 



dz, 



#o Jo (1 + zf v /l + fi ((l + ^) 3 -*+ 1 -l) 

and shall use this quantity to place bounds on the actual lookback time T(z) 
It is easy to obtain the leading term for f2 ~ 1: 

„ , s 2 



(52) 



(53) 



(54) 



3H (1 + w* 



-(3w*-l)/2 



} + o[n -i}. 



(55) 



The next sub-leading term again is trickier. We eventually obtain 
2 



T,, u (z) 



3H (1 + 



-(3iu,-l)/2 



} 



(1 + z )-3(t«.+l)/2 1 - (1 + z )-(9w*+5)/2 



(56) 



3(w* + 1) 9u>* + 5 

+ O ([Qo ~ I] 2 ) • 

Again, in any realistic situation (provided you accept the standard consensus cosmology) 
the uncertainties in w will completely dwarf any possible effect due to uncertainties in 
Qq. Exact integration and subsequent evaluation of the result for T Wt (z) can only be 
performed in terms of hypergeometric functions. Let us first be a little more careful 
about the use of the dummy variable in the integration and write 



T,,u (z) 



dz, 



Ho V^h Jo (1 + 5)2+1/2(3^+1) ^-(l-^^l + ^-^.+l) 

and then, (following the procedure of E]), apply the binomial theorem 



(57) 



[l - (1 - J^ 1 ) (1 + ^ ^ + T^ = £ ( ^) (-If (1+ 

n=0 ^ ' 



5) 



-(3iu*+l)n 



• (58) 



Now this particular binomial series will converge provided 



< 1. 



| (1 - Hp x ) (1 + 5)-(3«.+i) | 
That is, provided 

ll-^o 1 ! < (l + z) 3w * +1 . 



% Note that for z < one is actually calculating the "lookforward time" 
expands by an additional factor of l \ z \ in each direction. 



(59) 
(60) 

the time until the universe 
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< (1 + zf w * +l ] Vz e (0,z) or I E 0,0), 



1-O 



which is equivalent to 



n 



< min{l,(l + z) 3u, * +1 }. 



11 



(61) 



(62) 



In all cases, to ensure convergence at redshift zero, we must certainly have 

|l — ^ < 1, that is fi e(l/2,oo). (63) 

If z > Oand (3w*+l) ^ 0, (w„ ^ -1/3), or if z < and (3w*+l) 0, (w* < -1/3): 
Then (1 + z^ 3 ™* 4 " 1 ) > 1, and no additional limitation is imposed. 

If z > and (3™*+l) < 0, (w* < -1/3), or if z < Oand (3w*+l) > 0, (w* > -1/3): 
In this situation (1 + z)^ 3w * +1 ^ < 1, therefore we now obtain an additional limitation 
on z that is necessary to ensure convergence: 

— If z > 0, then we need 



z < 



Sk-1 



n 







If z < then we need 

1-Q 



z > 



n 







-1/|3id*+1| 



1/(3uj*+1) 



1 > 0. 



1 < 0. 



(64) 



(65) 



In view of equation (1221) these last conditions can also be interpreted as constraints 
on the Q parameter at the redshift one wishes to probe: 



< 1, 



that is Q Wt (z) G (1/2, oo). 



(66) 



Subject to this convergence condition we can integrate term by term, and obtain the 
convergent series 



T w . (z) 



1 M/2\ , ^ Jl-^T [!-(! + g )-(^+D"-3/2(^ + i)] 

Fov^^V " / (3w* + l)n + 3/2(w* + l) ' 1 ' 

As a practical matter, for many purposes this series representation may be enough, but 
we can tidy things up somewhat by first defining 

-1/2^ 



E 



n=0 



in which case 



T wA z ) 



n J (3% + l)n + 3/2K + l)' 
S Wt (1 - ^) - (1 + z)- 3/2K+1) S w 



(1 + Z )3».+l 



(68) 



(69) 
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Finally we Cciri recognize thctt (x) is itself a particular example of a hypergeometric 
function, and so we can write 

1 



3/2 (w* 



A3 


w;* + 1 


1 


9u>* + 5 


V2'2 


_3w* + 1 


' 2 


_3w„ + 1 



(70) 



Therefore 



T w , {z) 



3/2(^ + 1)^^ 



o 



A3 


u>* + 1 


1 


9io* + 5 


V2'2 


3w* + 1 


'2 


_3w* + 1 



; i - «o 1 



1 + ^-3/ 2 K + i) 



+ 1 



1 

'2 



9u>* + 5 
3u>* + 1 



1 ~ »o 1 

(1 + ^)3^*+l 



+ 1. ' 

As usual, T^z) can be used to bound T(z). Specifically, let w(z) lie in the bounded 
range w(z) G [u>_,u>+], then independent of fl ^ 1: 

• T w+ (^) <T{z) <T w _{z); (z>0), 

• T w _{z) <T{z) <T w+ {z); (-Kz<0). 

Here T w± (z) is given by the rather formidable equation (J7B . based on the use of 
hypergeometric functions. 



8. Special cases and consistency checks 

Useful special cases, and consistency checks we can perform on the formalism, include: 
Dust: For pure dust, w + = iu_ =0, we have simple exact results 
H dnst {z) = H (l + z)^l + n z. 

n (i + z) 



^dust(^) 



Vtn Z 



Pdust{z) = Po (1 + zf . 

, , . 2c r (vi + n oZ - 1)(V1 + ^ - 1 + «o) 



Tdust(^) 



+ 



F (l-Oo) 

^0 



(1+^ 

1 + 2 



(72) 
(73) 
(74) 
(75) 

(76) 



^o(l-fio) 3 / 2 
The only one of these equations for which the f2 
is the lookback time, for which 

2 f 1 



1 yr^rs v / r Tr s5J 



- dust, Cq=1 



1 



3H {' (1 + zf/ 2 

The classical hypergeometric series is given by 

(a) n {b) n x" 



1 limit is even remotely subtle 

(77) 



2-Fi (a, b; c; x) = ^ 



n=0 



,1 ' 



where (a)„ = a(a + l)(a + 2) (a + n — 1) is the rising factorial, or Pochhammcr symbol. This series is 
convergent for |a?| < 1. 
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Radiation: For pure radiation, w + = u>_ = 1/3, we have 



irradiation^) = H (l + z) y/ 1 + fi [( 1 + z) 2 - 1] . (78) 

_ ^ (l + z) 2 ^ 

"radiation^j — . . 2 ■ I 1 v ) 

1 - iZ + "0 (1 + z) 

Pradiation(-) = A> (1 + zf . (80) 



C y/l + ^oKl + ^-1]-! 



-^radiation \Z) 



1 



y/l+^0((l + z) 2 -l) 
1+2 



(82) 



The only one of these equations for which the f2o — > 1 limit is even remotely subtle 
is the lookback time, for which 

T ra diation,n =l(£) = jl — ^ + ^2 | ' ( 83 ) 

Cosmological constant: For pure cosmo logical constant w + = tw_ = —1. We then 
obtain (now as equalities rather than inequalities) what would for the NEC have 
been a set of bounds, such as those presented in [61 [71 El [H] . (That is, a nonzero 
cosmological constant is right on the verge of violating the NEC.) 

Furthermore, comparing with previous results in the literature: 

• For u>_ = — 1/3 one has 

H(z)>H (l + z). (84) 

This reproduces the SEC lower bound previously investigated in [51 [7J El EH]- 

• For u>_ = — 1 one has 

H(z) > H (1 + z)y/{\ + z)-i + in -l] [(l + ^-2-i], (85) 

whence 



H{z) > H y/1 + [Uo - 1] [1 - (1 + zf] = Ho v/n + [l-no] {l + zf. (86) 
This reproduces the NEC lower bound previously investigated in [21 [7J El EH] • 
For w + = +1 we have 

H(z) <H (1 + z)y/(l + zY + p - 1] [(1 + zY - 1], (87) 

that is 



H{z) <H (1 + z)y/l + Sl [{1 + zy-l}. 
This reproduces the DEC upper bound previously investigated in [61 [7J El EH]- 
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9. Conclusions 



In the absence of any detailed understanding of the precise nature of the cosmological 
equation of state p(p) it is useful to examine the question of just how much can be 
deduced with limited information. In this article we have worked in terms of the 
w-parameter w(z) = p/p, and we have used the idealized case of constant if* as a 
"template" for comparison purposes with more realistic w(z). Specifically: 

• For constant the explicit results for the density p w S z ) an d Hubble parameter 
H Wt (z) are well-known. The explicit result for the Q parameter Q Wt (z) is less well- 
known, and the explicit results we have obtained for the angular diameter distance 
dp w ^(z) and lookback time T Wt (z) appear to be both novel and significant. 

• More importantly we have seen that these idealized results for constant u>* can be 
used as the basis for general comparison results that bound the various features of 
the Hubble flow in the following sense: If we know that w(z) G [w;_,w;+] between 
redshift zero and redshift z, then for monotonically evolving generic cosmological 
quantities X(z) we have derived a number of rigorous bounds of the form 

X w± (z) <X(z) <X w ^(z), (89) 

where we have explicitly seen that the direction of the inequality depends both on 
the precise details of the evolution of X(z), and on the redshift range of interest. 

Finally we point out that all of our bounds have been explicitly calculated for all signs 
of the spatial curvature k G [—1,0,-1-1], that is for all flo (though we have restricted 
ourselves to the physically very plausible > 0). The bounds we have derived are thus 
both very general and very powerful. 
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